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Abstract. We obtain an asymptotic expansion for the tails of the random variable T = 
arg max„gK {A2 (u) — ) where A2 is the Airy2 process. Using the formula of Schehr [37] that 
connects the density function of T to the Hastings-McLeod solution of the second Painleve 
equation, we prove that as i -> 00, P(|T| > <) = Ce-t'^(*)t-"5/32(l + 0{t~^/^)), where 
ip[t) —t^ — 2i^/^ + 3t^/'', and the constant C is given exphcitly. 



1. Introduction and statement of the main result 

Directed polymers in a random medium (DPRM) were introduced by Huse and Henley [27] 
to describe domain walls in a ferromagnetic Ising model with random impurities, sometimes 
called a dirty ferromagnet. In the two dimensional Ising model, a typical domain wall is a 
lattice path in the plane, and DPRM in 1+1 dimensions is a statistical mechanical model 
whose states are such paths. For concreteness, we consider the square lattice N x N, and the 
graph with edges connecting nearest neighbors, so that the midpoints of the vertical edges 
have the Cartesian coordinates (i, j + 1/2) for some z, j G N, and the horizontal edges have 
the coordinates {i + 1/2, j) for some i,j G N. On each vertex, we place a random weight 
Sij independently from some distribution, and consider some set of allowable lattice paths F 
which originate at (0, 0) and always move up and to the right. We then define the random 
Gibbs measure on F as follows. For a path P G F, 



^i{P) = - exp 



where 



exp 



Per 



is the (random) partition function, and T > denotes temperature. For the set of allowable 
paths F, usually we consider 



F = F 



(m,n) 



{up-right paths ending at the point (m,n)} 
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in which case the model is said to have point-to-point geometry, or we consider 

r = r„ = {up-right paths of length n} , 

in which case the model is said to have point-to-line geometry. 

In the limit as T — > 0, the Gibbs measure becomes a delta function on the path with the 
greatest weight, and the randomness in the model comes entirely from the random weights 
Eij. This is known as directed last passage percolation. In the case of point-to-line directed 
last passage percolation with geometric or exponential weights on sites, it has been proven 
that the limiting fluctuations of both the energy of the maximizing path, and of the 
location of the endpoint of the polymer can be described in terms of the Airy2 process. 

The Airy2 process [M], which we denote A2{u), is a stationary process whose marginal 
distributions are the Gaussian unitary ensemble (GUE) Tracy- Widom distribution [ID], and 
is expected to be a universal process governing the limiting spatial fluctuations of random 
growth models in the Kardar-Parisi-Zhang (KPZ) universality class [29] in 1 + 1 dimension. 
This has been proven in the case of the polynuclear growth model [28l . Now let Ai and 
T be the random variables defined by 

Ai := max (A2(u) — u^) , 

u&R ^ ' ' ^ 

and 

T := argmax (^2(^) — ^^^) • 

Then M. describes the limiting fluctuations of the energy of the maximizing path in geometric 
point-to-line last passage percolation. The limiting fluctuations of the endpoint of the path 
are described by T. This fact was proved by Johansson [28] assuming that the maximum of 
A^iu) — is attained at a unique point almost surely. This assumption was later proved by 
Corwin and Hammond [12]. For DPRM at finite temperature, there are some recent results. 
For a continuum version and a semi-discrete version [33] of DPRM, as well as discrete DPRM 
with log-Gamma weights, the correct scaling exponents and a limit theorem for the free 
energy have been obtained at finite temperature [6l [H [TTl (311 [TOl [161 IM]- A limit theorem 
for the fluctuations of the endpoint of the polymer with point-to-line geometry, however, 
has only been proven in the case of geometric or exponential last passage percolation [28]. 
Nonetheless, T is expected to govern the fluctuations of the end-point of a polymer in DPRM 
at finite temperature as well, for a wide range of random weights Eij. For a review of the 
KPZ universality class in the physical literature, see [2S]- For a more recent review in the 
mathematics literature, see [2]. 

Remark: Although T has only been proven rigorously to govern endpoint fluctuations in 
the case of geometric last passage percolation, a formula for the endpoint fluctuations of a 
continuum version of DPRM was recently obtained nonrigorously by Dotsenko [20] • This 
formula is equivalent to the known formulas for T, given in fll.4p and (11. 101) below. The 
author of that paper was unable to prove this equivalence, and for the sake of completeness 
in the literature, we give a short proof in Appendix [O 

The distribution of AX is, up to rescaling by a constant, the same as the Tracy- Widom 
distribution for the Gaussian orthogonal ensemble (GOE) [H]. Specifically, we have 

P(7W <t)=Txifl'h), 
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where T\ is the Tracy- Widom GOE distribution function, defined below in (11. 5p and (I1.12p . 
This fact was proved by Johansson [28] by first proving a functional limit theorem for the 
convergence of the polynuclear growth (PNG) model to the Airyg process and using connec- 
tions between PNG and the longest increasing subsequence of a random permutation found 
by Baik and Rains [5]. A more direct, although nonrigorous, proof was given in [23] by 
analyzing the fluctuations of nonintersecting Brownian excursions. A rigorous direct proof 
based on the explicit determinantal formula for the Airy2 process was given in [17] , and the 
approach of [23] was made rigorous in [30]. Since the Tracy- Widom GOE distribution has 
been well studied over the past 15-20 years, a lot is known about the distribution of A^. In 
particular, the asymmetric tail behavior of T\ is 



5 l.|3/2 



O (s^^/^) ] , as s +00 

(1.1) 



^^^l^fl + 0f|.|-3/^)l, as. 



-oo 



62 



5C'(-1) 



where 

~ 211/48 ' 

and C(-) denotes the Riemann zeta- function. Similar formulas exist for the Tracy- Widom 
GUE and Gaussian symplectic ensemble (GSE) distribution functions, see [21 |3l HE]. Re- 
cently, some similar formulas have appeared for the general /3 Tracy- Widom distributions as 
well [m ESI [21]. 

Much less attention has been dedicated to the study of T ■ Exact expressions for the 
joint distribution of (A^,T) were obtained in two recent papers: in [32] by Moreno Flores, 
Quastel, and Remenik; and in [37] by Schehr. The formula of [32] involves the Airy function 
and the resolvent of an associated operator and is derived rigorously, while the formula of [37] 
involves a solution to the Lax pair for the Painleve II equation, and is derived nonrigorously. 
It was shown in [3] that these formulas are indeed the same, and therefore the formula of 
[37] is put on rigorous footing. 

Let us describe the two formulas. In order to do so, we first need to fix some notation 
which we will use throughout the paper. Let P(m, t) denote the joint density function of 
(A^,T). Let Ai(x) be the Airy function [7], and let be the integral operator acting on 
L^[0, oo) with kernel 

B,(x,y) = Ai(a: + y + s), (1.2) 
where s e M is a parameter. It is known that 1 — is invertible for any s G M. Let 

p,(x, y) = (1 - B,)-i(x, y), X, 2/ > . (1.3) 

Define, for t, m G M, 

'^{x\ t, m) = 26^*[tAi(t^ + m + x) + Ai'(t^ + m + x)]. 
Then the formula of [32] is 

/•oo roo 

P{m,t) = 2^'^7i{2'^'^m) I i ^{2^'^xu-t,m)p22/im{xi,X2)^{2^'^X2]t,m)dxidx2. (1.4) 

Jo Jq 
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In terms of the operator B,;, the GOE Tracy- Widom distribution function T\ equals 

J-i(s) = det(l - B,), (1.5) 

see [221 [36]. 

We now present the formula of [37]. To this end let g(s) be the particular solution of the 
second Painleve equation, 

g"(s) = sg(s) + 2g(s)^ 

satisfying 

g(s) ~ Ai(s) , as s — 7- +oo . 

This particular solution is known as the Hastings-McLeod solution [261 [23], ^-iid its uniqueness 
and global existence are well established. Consider now the following Lax-system (cf. [8]) 
associated to the Hastings-McLeod solution of the Painleve II equation, i.e. the system of 
linear differential equations for a two-dimensional vector $ = $(C, s), 

-^S*. (1.6) 

where the 2x2 matrices A = A{(, s) and B = B{(, s) are given by 

A(Cs)-( Ae + s + 2q' + 2q'\ 



B{C,s)=[ I ^ ) . (1.^ 



and 

-C -q 

The above system (11.61) is overdetermined, and the compatibility of the equations implies 
that q{s) solves the Painleve II equation. Now let $ = (^^) be the unique solution of (11.61) 
which satisfies the real asymptotics 

$i(C; s) = cos (^^C' + + O(r'), $2(C; s) = - sin (^^C' + + 0(0 , 

as C — ^ ±oo for s G M. Such a solution exists (see e.g. [H [191 [23]), and it further satisfies 
the property that $i(C; s) and $2(Ci ^) ^i-^e real for real ( and s, as well as the symmetry 
relations 

$i(-C; s) = <l>i(C; s), <I>2(-C; s) = -<I>2(C; s), C e C, s g M. 
Define for non-negative w and real s, 

POO 

h{s,w):= / C^2iC;s)e-^^'dC. (1.9) 



The formula of [37] for the joint density of Ai and T is 

P(m, t) = 4P(22/3m, , (1.10) 

where 

4 

P(s,w) = —J^i(s) / h{u,w)h(u,—w)du, (l-H) 
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and h{s, —w), w > is understood as the analytic continuation of w) to the negative real 
axis. In terms of the Painleve function q{s), the Tracy- Widom GOE distribution function 
J-'i can be written as HT] 



2 roo 1*00 roo 

J'i(s)=exp -- J q{x)dx - - J J q{xf dx 
The marginal density function for T equals 

/oo 
P{m,t) dm, 
■oo 



dt 



;i.i2) 



P{t) ■= 2^''^P{2^'h) , (1.13) 



P{w) = / P{s,w)ds. (1.14) 



and using (11.101) . we have 
where we introduced 



Our main result in the present paper is an asymptotic expansion of the density function P(t) 
for large t. 

Theorem 1.1. As t ^ oo, the marginal density function P{t) satisfies 

P(t) = re-l^(t)t-^y32 + + 0(^-3/2^ 

which extends to a full asymptotic series in powers of t~^^^. The function ipit) is 

ifit) =t'- 2t'/^ + 3t3/^ 

and the constant r is given by 

r = 2-29/6g5/4giC'(-l)7r3/2^ 

with ({z) denoting the Riemann zeta- function. 

Integrating the above density, Theorem 11.11 implies the following corollary. 
Corollary 1.2. As t oo, 

P(|r| >t) = Ce-t^Wt-i^5/=^2 A + + 0(r 3/2)^1 . (1.15) 



where C = r/2. 



Remark: The leading decay order of e~'^^^ in (11.151) was first predicted in the physics 
literature by Halpin-Healy and Zhang [25], compare also [31]. The first rigorous confirmation 
of this rate of decay appeared in the paper [15] of Corwin and Hammond, in which they 
give e~'^*^ as an upper bound on P(|T| > t), although they do not give the value of the 
constant c. In the paper [37] of Schehr, the author uses the formula (11.101) to find that the 
leading coefficient is c = 4/3, although he did not employ complete and rigorous estimates. 
Recently Quastel and Remenik [35] rigorously obtained explicit bounds on c, and remarked 
that they believe that the correct rate of decay is in fact 4/3. Corollary 11.21 confirms this 
rigorously, and in addition gives subleading terms and constants. In principle, all terms in 
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the asymptotic expansion are computable by the methods of this paper, but calculations 
become more involved. 

The setup for the remainder of the paper is a follows. We will prove Theorem 11.11 using 
(11. lip . (I1.13P and (I1.14p . In Section |2] we evaluate h{s,w) as w — oo by analyzing the 
system (II. 6p close to the origin. In Section [3] we prepare h{s, —w) for asymptotic analysis 
as w — oo using the global identity from [4] which allows us to express $2(Cj ^) terms of 
the resolvent of a Hankel operator on L^[0, oo) whose kernel is constructed out of the Airy 
function. Then in Section |4] we split P{w) into two parts, one which can be estimated using 
asymptotics (II. ip of the Tracy- Widom distribution function J^i(s), and another which we 
evaluate asymptotically by Laplace's method. Theorem 11.11 then follows via (I1.13p . and the 
proof of Corollary 11.21 will be given in Section O 



2. Expansion of h{s,w) as w ^ +oo 
With the change of variables A = CV^w, (II. 9p becomes 



h{s,w) = ^[ A<l>2 ) e-^rfA. (2.1) 



We thus see that we need the Taylor expansion of $2(Cj"S) at ( = 0. To this end let us 
consider (11.60 as a system for a 2 x 2 matrix- valued function 



9* 

9C 



(E^n(.)c)^, ^ = (E^n(^)cV, ncs) = mcs),Hc,s)), (2.2) 

^ n=0 ^ ^ ^ n=0 ' 



where 



, t;(s)+V(s)\ ,_Ag(s) \ , 4 

- t v{s) J ' - -4g(s) J ' ^2 - Q 



-q{s)J ' ' \-l 

are determined from (ll.7p and (II. 8p . and 

v{s) := s + 2q{sf -2q'{s) . (2.3) 

The matrix function can alternatively be defined as the solution of a certain oscillatory 
Riemann-Hilbert problem, see Appendix \M Such Riemann-Hilbert formulation appears in 
one form or another throughout the literature on the integrable structure of Painleve II, and 
is the basis for many results (cf. |23]). In particular, we will use the facts that 

^(-C; s) = as^iC; s)as, (eC, seR, (2.4) 

and 

^(0;s) = e~"3^=°°«W'^*, (2.5) 

which are justified in Appendix |Al 

Since the coefficient function A{(, s) is analytic at the origin, there exists a solution \E'(C, s) 
of equation (12. 2p which is holomorphic in some neighborhood of the origin and uniquely 
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determined by its value \l/(0,s) = ^I/qI-s) (cf. |12])- The Taylor expansion of \i/ about C = 



IS 



^{C,s) = J2'^^{s)C, \C\<p, p>0. (2.6) 



n=0 



Inserting (12.61) into (12.21) and comparing coefficients yields for n G Z>o &st 

{n + = Ao^nis) + Ai^„_i(s) + A2^n-2{S), ^_l(s) = ^_2(s) = 0, (2.7) 

and secondly 

^ = g(s)(T3^'„(s) +2a2^„_i(s), (^3=Iq _i 1.^2=1^ q1- 
The latter differential recursion implies 

for some 2x2 invertible matrix valued constant Cq. In fact, from (12. 5 p we see that Cq = I. 
In order to obtain the rest of the Taylor coefficients, we use the recursion (12. 7p to write 
explicitly in terms of the matrices Aq, Ai, A2. The ffist two coefficients are 

^i(s) = Aoe-'s/r^W'^^, ^2(5) = ^(^0 + Ai)e-"3^''°°^(^)'^^ (2.8) 

At this point we would like to connect the coefficients "^nis) to the required values of $i(C; s) 
and $2(C; ■5). The aformentioned global symmetry relation (12. 4 p implies 

^ni-s) = (-l)V3^'„(s)a3, n e Z>o, 

i.e. all coefficients '^2n{s) are diagonal whereas all coefficients \l/2n+i(s) are off-diagonal. We 
therefore have 

00 00 

?i=0 n=0 

(2.9) 

From ([22D (see explicitly ([2l])), 

=^n(s,g(s),g'(.))e-r^W<^-, (2.10) 



with polynomials Pn in three variables. Integrating (12. ip . we obtain the following Proposi- 
tion. 

Proposition 2.1. As w ^ +oo, the function h{s,w) satisfies the asymptotic expansion 



Ms,^) = £^e-r^W'^-^Q„(,,g(,),g'(,))^-, (2.11) 

n=0 



g„(.,g(.),g'(.)) = er .(^)-v,(-),(.)(^^, 

is a polynomial in 3 variables. This expansion is uniform in s on compact subsets of the real 
line. The first two terms in the expansion (12. lip are 

Qo{s,q{s),q'{s)) = -v{s) , Qi{s,q{s),q'{s)) = ^-i^ + v{s){v{s)q'{s) + q{s)) -2 , (2.12) 
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where v{s) is defined in i \2.3L Moreover, as s ^ —oo, the polynomial Qn{s,q{s),q'{s)) 
satisfies 

Q„(s,g(s),g'(.)) = 0((-s)("-i)/2). (2.13) 
Proof. Replacing $2 in (12 .ip by its power series (12. 9 p gives 

00 ,^,(21) 



^^"'^^ (2^)3/2 (2^)n I 

00 



n=0 



-W 



w 

(21) 



Using the structure (I2.10p of ^I^2n+i('5)5 can thus write, as w — )■ 00, 

h{s,w) = ^e-i'^°°«(-)'^-^g„(.,g(.),g'(s)V- 

n=0 

For the asymptotics (12.131) . we use [19^1 . Namely as s — )■ —00, 



, -s / 1 73 10657 ^, _i2, 
^^^) = ^/^(^ + 8i^-128i^+ 1024i^ + ^('^ 

1 / 5 803 181169 ^, _i2^ 
^^^) = -WV^ 1^-8^^+1287 -1024;^ + ^^'' ' 



(2.14) 



from which it follows 



1 (-s)-3/2 5 9(_3)-9/2 ^ 1323(-s)-i5/2 



2V2 8s3 4^2 128s6 32 

181169 



1024s9 



+ 0((-s)"2i/2) 



From these we find that, as s — >■ —00, 



(s)-( 0((-s)("-l)/2)\ ..^(,),, 

VE^2n+l(s) - ^^((_^y„_l)/2) Q je s 

which is easily checked by induction using the recursion (12. 7p . □ 
Let us now consider h{s,w) as w — )■ —00. 



^We would like to point out the following small detail. In [2] the coefficient of 0(s ^) is given as 



10219 
1024 ' 

whereas in [H] and gUj it is 
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3. Airy formula for h{s,w) and discussion 

In order to analyze h{s, w) for large negative values of w it is convenient to use a formula 
for $2 which was obtained in [1]. Let us review the notation of j^. As in (11 ■2p . let be 
the integral operator acting on L^[0, oo) with kernel 

Bs{x,y) = Ai{x + y + s). 

Let As := be the Airy operator acting on L^[0, oo), which has kernel 

POD 

As{x, y)= Ai(x + s + OAi(y + s + ^)di 
Jo 

_ Ai(x + s)Ai'{y + s) - Ai'{x + s)Ai{y + s) 
x-y 

Define the functions Q and R as 

g := (1 - As)-'Bs6o , R:={1- A^y^AJo , (3.1) 

where 6o is the Dirac delta function at zero. As in [3] we use the convention that the Dirac 
delta function satisfies J^Q^^6o{x)f{x)dx = /(O) for functions / which are right-continuous 
at 0. Introduce also the function 

e(x) := - sin ( + (s + 2x)C 



Proposition 2.1 of |1] is 

$2(C, s) = e(o) + (9, i? + g), (3.2) 

where the functions Q and R are defined in (13. ip . and (■, ■) is the inner product on L^[0, oo). 
Notice that the function R + Q can be written in terms of the resolvent kernel of the operator 
B5. Indeed, let 

R, :=(1-B,)^-1, 

be the resolvent of the operator B^. It is an integral operator as well, and let us use Kg to 
denote the function Rs^o- We then have 

R + Q = {1- Bl)~\Bs + Bl)6o = (1 - B,)-iB,5o = (1 - Bs)-\l - (1 - B,))5o 

= {{l-B,r'-l)6o = Ks, ^^'^^ 

and (13. 2p becomes 

$2(C,s) = e(O) + (0,A',). (3.4) 
Substituting (13.40 into (II. 9p gives 



with 



his, w)= - j\ sin + e-'^^'dC + C (0, Ks) e-^^'dC 
= a{s,w) + J a{s + 2x,w)(^{l -Bs)~^TsAi^{x)dx , 

a{y, w) = - j\ sin + e'^'^'dC, 2/ G K, 



(3.5) 



(3.6) 



and the translation operator 

{Tj){x) = f{x + s). 
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When w lies off of the non-negative real axis (13. 6 p diverges. In order to obtain the analytical 
continuation of a{-,w) to the complex w-plane, set ( = —iz. Then 

a(y, w) = - z sinh -^^ - yz e'^^'dz = - / 2e3"'+"'"'-^"rfz, 

and by Cauchy's theorem we can deform the contour of integration such that 

a{y,w) = ^ [ ze-^''+'"''-y'dz, (3.7) 

where £ is any contour that starts at the point infinity with argument arg z = (see 
Figured]) and ends at the conjugated point. We make the following observations. First (13. 7p 




Figure 1. The integration contour C chosen in (13. 7p . The dashed lines are 
the rays arg z = ±| and arg z = ±|. 

converges absolutely and uniformly in any compact w domain, hence (13. 7p determines the 
analytic continuation of a(-, w) to the entire real line and it allows us to consider (13.51) now 
for all w G M. Secondly using the definition 

Ai(z) = — [ e^''~''dt, zeC, 
2m Jc 

we obtain (see |1] for a similar computation), 



aiy,w) = ^e^+^ ( Ai'(2-^/^y + 2-«/V) + ;^Ai(2-^/^y + 2~'/W)], w,y eR. 



w 



In particular, notice that, for y ^ —w"^ /A we have, as — oo. 



a{y, -w) 



g 24 4 



1/4 
n 



E 

n,m=0 

-1/4 



4y 



n 



exp 



3 
2 

2 
3 



3m 
2 



(3.8) 
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with (cf. [7]) 

c-d-1 , _^X3n+j^ 6n + l 

We now go back to (13.51) and investigate the kernel Ks{x). We consider it as a function of 
both X and s, so let us write 

K{x, s) = Ks{x) = ((1 - B,)-iT,Ai) (x), x> 0, 

which solves the Fredholm integral equation 

/•CO 

Ai{x + s) = K{x,s) - / Ai{x + y + s)K{y,s)dy. (3.9) 

Jo 

Recall at this point that the Hankel operator : L^[0, oo) — ?■ L^[0, oo) is a bounded 
operator and 1 — is invertible for any fixed s G M. Thus K[x, ■) G L^[0, oo). In fact, as a 
consequence of (13. 9p . we have K{x, ■) G C°°[0, oo) fl L^[0, oo) for fixed s G M and 



K(x,s) = Ai(x + s)(l + o(l)), 



X — > oo, 



uniformly on any compact subset of{s : — oo < s < oo}. For regularity of K{s,-) as a 
function in s, we get via Holder's inequality and (13. 9p . 

K{s,-) G C°°(M). 

To obtain a statement on the large positive and negative s behavior of K{s, ■), we can view 
Bs equivalently as an operator B acting on oo) with kernel 

B{x,y) = Ai{x + y). 

A small norm argument for B then shows that K{s, ■) is approaching zero exponentially 
fast as s — 7- +00 for any fixed x G M. On the other hand from (13.91) . K{s, ■) can grow at 
most power-like as s — )■ — oo, again for any fixed x G M. For our purposes we are especially 
interested in the behavior of K{0, s) as s — oo. By (13. 3p . K{0, s) is 



and thus (see pjli]; 



Hence (cf. [19], [lO]) 



K{0,s) = R{0) + Q{0), 

q^{t)dt + q{s). 



4. Large w expansion of P{w) 

We have now gathered enough information to compute the asymptotics of P{w) as w — )■ oo. 
Let us write P{w) as 

P{w) = Pi{w) + P2{w), 



12 

where 
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h{u, w)a{u, —w)du 



ds . 



h{u,w)< / a{u + 2x, —w)K{x,u)dx > du 



J — oo 

/oo 
J-i(s) 
-oo 

We first focus on Pi{w). 

4.1. Expansion of Pi{w). Using (12.111) . we liave as w — > oo, 

,3/2 



ds. 



Pi{w) 



-e 24 



2l0/3y;3/2 

where 

Gn{s, w) = Qn{s, q{s), q' (s)) (Ai'(2-2/3s + 2-'/V) 
Let us write this as 

^3/2 



/oo ^ 1 /"OO 

■oo n=0 ""^ 



W 

2^73 



Ai(2-2/3s + 2-8/3y;2^y (4.1) 



AH 



-e 24 



2l0/3y;3/2 ^ 

n=0 



'Hn{s,u) duds , 



(4.2) 



with 



7^4s,m) = J-i(s)e-^-i''^"''(^)'^^G„(n,ti;); {(s,n) G 

At this point recall the asymptotics from 0, as s — oo, 

Ti P. 24 



-oo 



< s < u < oo}. 



Tie 24 3^ 



11/16 



55|s|-3 l0675|s|~9/2 
24^2 2304 165888^2 



+ 



3970225|s|- 
31850496 



as well as for u — )■ — oo. 



-3/2 



\u\ 



+ 0{\sr'/')), n 



2629|m|-9/2 105131 



211/48 ' 



(4.3) 



u\ 



I2V2 576 20736\/2 1990656 



Since Gnl^, u') has at most polynomial growth as s — — 00 and is exponentially decaying as 
s — 7- 00, we see that 'H„(s, u) is indeed integrable over (here we used that J-'i{s)e~ ■l^^ g{x)dx 
is bounded as s, u — )■ -foo). Thus we can apply Fubini's theorem in (14. 2 p and obtain 



nn{s,u)duds= / F{u)e-'^-^^'i^''^'^^Gn{u,w)d 



(4.4) 



where 



For large negative u, we can integrate (14. 3 p and obtain the following proposition. 



F(u) := / Ti(s)ds. 
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Proposition 4.1. As u ^ — oo, 

, 8ri _Hi_k^/ 97^2 19337 24666605^2 

F(y] = 1 p 24 3^2 I 1 1 

^ ^ |n|33/i6 \^ 48|u|3/2 2304|n|3 331776|n|9/2 

1358238769 ,„i5/2n\ 

H r^ + 0(\u\ ^^/^) , (4.5) 

31850496|m|6 ^' ' ^J^ ^ ' 

and Ti is given in (14.31) . 

Proof. We start with (14.31) and obtain for \u\ sufficiently large after a change of variables 



3 3/2 

~ e"^"Wf / 55x-3 10675X-9/2 



^^""^ 'lui x^l^^ V 24^2 ' 2304 165888^2 



3970225X 
31850496 



-6 



+ 0(a:-^^/2))rfa;. (4.6) 



Put ^ = + or equivalently, 



so that as t — )■ oo, 



x = »--(-l + v^rT3t)'^', (4.7) 



X 



^ 3^3 27t 243^3 2187^2 ^ > 



and 



24 V3/ t-i/2 4 35^3/2 80 \ 

3^2/3 V 3^3 27t 486v^ 2187^2 ^ >) 
Substituting (14. 7p into (14. 6p . we obtain 

124^/1^ , 35ri/2 8229 1278127r3/2 



F{u) = ^ r e-V^i/i*^fl 



96^3 55296t 15925248^3 

661165345 ^, , 

which can be readily integrated by parts to produce (14. 5p . □ 
Notice in particular that by (14. 3 p and (14. 5p . as m — t- — oo, 

JlH ^!f!A,^,J3__ 97V2; _ ^ . i5/2nA 

F(n) 8 V |^^|3/2 ^ 2|n|3 4|u|9/2 16|n|6 ^ ' V' 



We split the integral in (14. 4p in the following way. For some a > 0, which will be specified 
more precisely in the following, let 

Hn{s,u)duds = Ii^n + h,n, (4.8) 
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where 

-oo 

/oo 

For choose w > sufficiently large and replace F{u) by its asymptotics (14.51) : 

■oo ^ ^ 

e-[-'^^--^^-^'^''^''Vt~^V^^G„(-t,w)(^l + 0{^^^^^ (4.9) 

Since (see [2]) 

j ^ q{x)dx = ^t^/^ + ^log2 - J (^q{x) - )J^^dx, 
we can use (I2.14p and obtain 

1^ qi^^)dx = ^£3/2 + ^ log 2 + + O , t oo. 

To evaluate (14. 9 p asymptotically, we use these asymptotics to obtain 

/"OO / ,3 ,3/2 ,\ 

/i,„ = 4v^ri / e-^-+— -^)r=^='/i6G„(-t,ti;)(l + 0(r3/2))rft. 

In order to ensure convergence of the integral, we need to impose 

24 _n/n 6w 
1 H i=t ^ > 0, as w oo, 

which can be guaranteed for any a > |. Let us use the change of variables 

u = — -\ ^ -, t>0. 

24 V2 4 ' 

As M — )■ CXD 

t = 24^/3^^/3 A _ + ^1^^-2/3 + o(n-i) Y (4.10) 



3 12 

and we notice that the only critical point of the function u = u{t) as — )■ oo is given by 

to = v^(l + 0(i/;-3/4)), 

which by the choice a > | lies not in the domain of integration. Moreover, 



3 V 6 12 V ; ; 



and we conclude 



h,n = ^V2n24-''/'' r e-u-''/''Gr.{-t{u),w)(l + 0{u-'/'))du. (4.11) 



24 ^ ^ 
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Since Gn{—t{u),w) in the integral above grows at most like a polynomial as w — oo, we 
have the estimate 



\h,n\<ce 24 V2 + 4 \R^(^w°')\, as w cx), (4.12) 

for some polynomial i?„ and c G M. 
For /2,„, 

/oo 

we choose more specific bounds on a, namely 1 < a < 2. Then, as w — ?■ oo, 

2-2/3^ + +00, 
and we can use the expansion (13.81) for a{u, — w). This gives 



Gn{u,w) = -2-2/3^/^Q„(n,g(«),g'(n))^(u,^)e-i(2-''''^+^-'''-')''^ ^ ^ oo. 



where n(u, w) can be read from (13.81) 

°° / /( \ " / O \ ™ 3m 



1 

= ^ + 3^ + 2^ + ^^^"')' ^^-^^^ 

In general, the coefficient of w~'' in the latter series is a polynomial in u of degree at most 
k/2. Let us now write 



IT 

where 

/oo 
-w°' 

and 

1 ^ 1 .,,2 / ,„a3/2 



iJ(u) = iJ(n|tf;) = lnF(u)H h— / (/(x)dx H (1 + — | ,nG 

w 4 w 

The derivative of the function Hiu) is 



i^'(n) = -4^^ + 4-^ + -,/l+'" 



wF{u) 4 u; 4 V 

As — oo, the only zero uq of the function H'{u) will lie in a neighborhood of n = — oo, so 
we can use Proposition 14. ![ (14. 3 p and (12.141) to determine it. The solution is 

/ 3 65 3 \ 

«o = -2V3;(l-^- 5^-5^ + u.^oo, (4.14) 
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and the integral Ln{w) can be evaluated as w — oo by Laplace's method. For the convenience 
of the reader, we have included a short review of Laplace's method in Appendix [Bl Notice 
that as w — >■ oo, 

_iimo_u£^ , 41^")^/^ 3/2_3 3/4_97 / 21 , „Q/2\A 

-r°°q(x)dx -4^3/4/ 35 , -3/2\\ 



V2 V 12^3/4 

which implies 

e-^Kk) = ^'^''"l''|:'^"'^' e-#+l-^^^-^-^^^ f 1 + ^^-3/4 + 0(^-3/2) Y (4.15) 
^dd/d^ y 2 / 

Also, 

M"W)-"^ = ^(l+3|7lH-0(»-/^)). (4.16) 
Combining fl4.15p and fl4.16p . Laplace's method gives 

229/54y^^iC'(-l),5/4 _^^,^3/2_2^3/4^ . , , V ^^ 



It follows that 



913/24 iC'(-l)g5/4 3 

r ^ C _iu 1 2 3/2 _ 2^,3/4 / / \ // \\ 

/2,n = e 24 +3^^ Qn{uo,qiuo),q{uo)} 



w 



25/32 



xn(no,u;)(^l + ^^ + 0(ty ^ 00. 

From f l4.13p and the property that the coefficient of w~'^ in the asymptotic series for Il{u, w) 
is a polynomial in u of degree at most k/2, we find 

U{uo,w) = l + 0{w-'^), w ^ 00. (4.17) 

Plugging into 04.21] and fl4.8p . we obtain 



^3/2 
21073^3/2 ' 



n=0 

1 A'/ 



00 

X 

n=0 



73/32267/24 I ' §^3/4 

oo 

Xg„(uo,g(uo),g'(uo))u;-", (4.18) 
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where we have used the estimate (14.121) which shows that Ii^n does not contribute to the 
leading order asymptotics when a > 1. Next by Proposition 12. ![ as w — oo, 



n=0 

and from f l212|l . 

2wV4 V 8^3/4 



If 5 



as — 7- OO, 



and 



Thus 



Qi (no, g(no), g'(no)) = + O {nj-^/^) . 



E gK), q'{uo))w- = ( 1 + ^ + 0(n;-='/^) ] , w^oo. (4.19) 

n=0 ^ ^ 

Combining this with (I4.18P we find, as — > oo, 

Pi(n;) = Kc-'^H-'^'-'-'^'w-''/'' (l + ^ + 0(n;-^/^)) , (4.20) 

with 

K = 2-91/24^3/2giC'(-l)e5/4. (4.21) 

We now turn our attention to P2. 

4.2. Expansion of P2{w). Following the philosophy of the previous subsection, write 

3 °° 1 f f f 
^2H = 2W3^3/2 ^"^$Z^ L'Hn{s,U,x)dxduds, 

n=o J J Jn 



with 

it?(ii + 2ac) TOO 



G,(M,M;,a;) = g„(M,g(M),g(M)) Ai — ^ + _ Ai — ^ + \ K{x,u) 



+ 2x \ w'^ fu + 

22/3 28/3 ] ~ 2^/3 I 2^/3 ' 2*^/3 

and the domain of integration 

Vt = {{s,u,x) e : -00 < s < u < 00, X > 0}. 

Recalling the behavior of K{x, u) as we approach the boundary of fi, we can again apply 
Fubini's theorem to obtain 

/oo ^ 
F{u)e-'^-^^ 5(*)'^*G„(m, w)du , 
-00 
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where (compare (14. ip ) 

Gn{u,w) - 



Gn{u, w, x)K{x, u)dx 
Qn{u,q{u),q'{u)) 



wx 

e 2 



Ai' 



f u + 2x w 



+ 



w 



:Ai 



u + 2x w"^ 



K{x, u)dx. 

(4.22) 




22/3 ' 28/3 j 2^3 \ 2^/'^ 2'^/'^ 

Introducing 1 < a < 2, we again split the integral as 

Tinis, U, x)dx du ds = /i,„ + l2,n 

' POO 

J-w°' 

For we follow the same steps as for and conclude (see (14. lip ) 

4. = |v^n24-iV48 r e--u-''/''G^{-t{u),w))(l + 0{u-'/'))du, 

24 V2 4 

with t = t{u) given in fl4.10p . Also here, Gn{u,w) grows at most like a polynomial as 
u —7- — oo, hence we have the estimate 



as w — )■ cxo, 



for some polynomial Rn- 

For J2,n we can replace the Airy functions in (14.220 with their large argument asymptotics, 
giving 

iw 



G^{u,w) = -2-^/' 



with 



J{u,x) = J{u,x\w) = 1 + ^(1 + 



^-u,j(u,x\w)Yi^^ + 2x, w)K{x, u)dx , 
A{u + 2x)Y^^ 



Plugging this into the integral l^^n gives 



'2,n 



-2-/^r 



-tu(H(ti|ui)+ J(u,x|ui)) 



D 



Qn [u, q{u) , q'{u))Il{u + 2x, w)K{x, u)dx du . 



where the domain of integration D equals 

D = {{u,x) e : -w" < M < oo, X > 0}, 

where a G (1, 2) is fixed, and H = H{u\w) is given by 

1 



u 1 

H{u\w) = In F{u) + - + - 

w ' ^ 2 



q{x) dx. 
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Since 



d " 

— {H{u\w) + J{u,x\w\) 
ou ^ 

d - 

— {H{u\w) + J{u,x\w\) 
ox ^ 



wF{u) 



1 _ qiu) 
4 w 



4(m + 2x) 



1 1 / 4(m + 2x) 



we see that H{u\w) + J{u^x\w) has no critical point in the interior of D. However its 
partial derivative with respect to u vanishes in a neigbhorhood of (— oo,x). More precisely, 
it vanishes for u = uo{x), where 



Uo[X) 



2v^ 1 



65 



2^3/4 32^3/2 



X 

+ — 



3 3x 
+ 



as w — oo. 



Notice that no(0) = uq as given in f l4.14p . Laplace's method now indicates that we need to 
expand the integral l2,n in a neighborhood of the point (mo(0),0). First, 



H{u\w) + J{u,x\w) = H{uq\w) + Jx{uo,0\w)x 

+ ^{u - Uo,x)M{uo,0\w) 



U — Uq 
X 



+ 0((^i-Mo)V), 



Jxx{Uo,0) = ^( 1 + ^ 



2v^ 



1 



1 



0/ 4^3/4 10 Oy 2^3/2 V 4 8 



-7 4 



9 



-1/2 



which is valid in a neighborhood of (uo(0),0), with M{u,x\w) denoting the Hessian matrix 
of H{u\w) + J{u,x\w). Since the mixed partial J^u differs from J^x by a factor of 1/2, we 
have 

M{uo,0\w) 
and therefore, 

M{uo,0\w) = 



1 



64^9/4 10 



as w oo. 



Here we have used that (compare fl4.16p ) 



H"iuo) 



1 



2v^ 



7 



9 



4^3/4 2m;3/2 64^9/4 



as w — 7- oo. 



Moving on with Laplace's method, we have 
Iw 



'2,n 



-2-2/^/^ 

TT 



D 



w{H{u,w)+j[uM^))q^ (m, q{u),q'{u))ll{u + 2x, w)K{x, u)dxdx 



-2-2/3^^^g-«,//KI-)g^^^^(^^)^^'(^^))n(Mo,w)K(0,Mo) 



TT 



-wJx(uofl\w)x^^ 
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The integrals in the above expression can be expanded for large w as 



and 



Combining with f l4.15p . we obtain 



2l3/24giC'(-l)g5/4 ^3 



2,n = e 24 +3"' Qn{uo,q{uo),q{uo)) 



w 



57/32 



xU{uo, w)K{0, uo) (^1 + ^ + O {w-'/')y 



Applying the same arguments as in the computation of Pi{w), we get from f l4.17p and f l4.19p 

-3/2 - °° 



n=0 

_Hi+2^3/2_2^3/4 



= 210/3^3/2 ^"^^ Z^^Ul." +-^2,nl 

-^/32jaO,«o) (l + ^ + , (4.23) 

with K, given in (14.21 p . At this point recall (I3.10p . which gives 

K{0,uo) =w(^l- -^ + 0{w~^/^)^, asw->cx). 

Inserting this asymptotics into (I4.23p . we find 

P,{w) = Ke-^+i'^'''-2-''"u;-«i/=^2 ^1 ^ _i3_ ^ ^(^^_3/2)^ _ (4 24) 

4.3. Asymptotic expansion of P{w). Since P{w) = Pi{w) + P2{w), we add (I4.20p and 
(H^ . which leads us to 

with the constant k 

^ = 2-9V2V/2g|c'(-i)e5/4_ 
Plugging into (I1.13p . we get Theorem 11.11 

5. Proof of Corollary 11.21 

Since P{-,t) is an even function, see (ll.4p or (ll.lUp and (II. lip , we use the symmetry of 
the density function P{t). As t — )■ oo, 

/oo 
P{s)ds 
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Let US use the change of variables 

U = ^{t) 2^3/2 + 3^3/4^ 

hence as n — )■ oo 



t 

and we obtain 



P(|r| > t) = f e-K"i^^/9^(^l + ^ + o(«-i/2)^dw. 

The latter integral can be readily integrated by parts and we obtain (11.151) . 

Appendix A. Riemann-Hilbert problem for $ 

As mentioned in Section [21 the matrix function \1/(C; s) described in that section can be 
defined as the solution to a Riemann-Hilbert problem (RHP). The exact formulation of this 
RHP is as follows. Define the rays Pj for j = 1, 2, 3, 4 as 

Pi = {C : argC = 7r/6} , P^ = {( : argC = 57r/6} , 

P3 = {C : argC = -57r/6} , P4 = {( : argC = -7r/6} , 

and give them the orientation from left to right, as shown in Figure [2l The plus- (resp. 
minus-) side of the contour is then the side which is lies to the left (resp. right) side of the 
contour when facing in the direction of the orientation, as shown in Figure [21 The RHP 
consists in finding the piecewise analytic 2x2 matrix function \1/(C; s) = {^ij{C] s)).j_^ such 
that 

(1) *(C; s) is analytic for C G C \ 

(2) For C G U^^^Pj, \I^(C; s) takes limiting values from the plus- and minus-sides, and 
these limiting values are related by 

vl/+(C;s) = vi>_(C;s)j*(C), 

where 



MO 



2 
1 

2 3 



(3) As C ^ 00 

^(C-s)- (l + 0(C-'))( cos(|C' + <) sin(|C=^ + <)\ 

vi/(c,sj-(i + u(c (^_gin(|c3 + sc) cos (|c=^ + • 

(4) \1/(C, s) is bounded for ( close to the origin. 

This RHP is uniquely solvable [9], and the vector $((^, s) defined in (II. 6p is the first column 
of ^. That is, 

$i(C;s) = ^An(C;5), $2(C;s) = v^2i(C;5)- 
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Since the jump matrices satisfy the symmetry relation 

A 0\ 

we also have 

*(-C;5) = ^3^(C;s)a3, (eC, seR, 

which is 

To justify (12.51) . we appeal to a result from the paper [S] of Baik and Rains. In that paper 
the authors consider a different RHP for a different Lax pair for the Painleve II equation. 
Our RHP can be transformed to theirs by a series of explicit transformations. After doing 
so, equation (12.51) is equivalent to Lemma 2.1 (iii) in [5]. 



Appendix B. A review of Laplace's method 



Let us give a quick review of Laplace's method. Recall that Laplace's method can be 
applied to the integral 

when the function H{u) has a unique minimum at Uq G {a,b), where a and/or b may be 
infinite and /(mq) 0. If we assume that the functions H and / are analytic, then Laplace's 
method gives that, as w oo. 



X 



5 fH'%uo] 



2wH"{uo) 



12 \H"{uo) 



-wHjuo) /(%)V^ 

AH"iuo) fiuo) 3/(mo) H"iuo)fiuo) i 



0(w' 



Another situation might arise if H{u) has no critical point in the interior of (a, b) but a 
global minimum is achieved on the boundary, say at mq = a with H'{a) > and /(a) ^ 0. 
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In this case Laplace's method gives as w -> oo, 



-wH{u) 



f{u)du 



wH'ia) 



1 

1 + - 



/'(«) 



H"{a] 



1 



/"(«) 



w \f{a)H'{a) {H'{a)Y 
3H"{a)f'{a) H"'{a 



w'\f{a){H'{a)y f{a){H'{a)r {H'{a)y 



Both of these cases can be generahzed to iterated integrals with appropriate assumptions on 
H and /. 

Appendix C. The formula of Dotsenko 

Using the so-called replica Bethe ansatz, Dotsenko recently derived a formula for the 
distribution of the fluctuations of the endpoint of a continuum directed polymer in a random 
medium. In this appendix we show that his formula coincides with the known formulas of 
[52] and [H7] for T, given in (ll.4p and fll.lOp . respectively. In his article [20], Dotsenko derives 
an expression for the distribution function W{x) for the fluctuations of a directed polymer's 
endpoint. The conjectured relation between W{x) and the density function for T is 

fOO POO -I 



W{x) 



P{w)dw, or equivalently W{x) 



4 



Pi2~'^''^m,2-^'h)dmdt, 



(C.l) 

where P{w) and P{m,t) are defined in fll.l4p and (11.41) . respectively. The formula of Dot- 
senko is (see equations (6), (7) and (8) in [20] ) 



W{x) 



POO 


POO P 






/ — oo 


Jo Jo 



Ps{uj,uj')^:^>^{s,x) dudu' 



ds 



(C.2) 



where J^i is the Tracy- Widom GOE distribution function, and and ps are as defined in 
(ll.2p and (II. 3p . respectively. The function ^^^i is given as the integral 



^uuj'{s,x) 



d_ 

duo 



d_ 
'du' 



-]<^{uj + - + y- x)^{uj' + - + y; -x) 



+ 



d d 
du du' 



^(wH y;x)'^{u' -\ \-y;-x) 



dy-. 



with 



'^(u:x) 



1 

27ri 



exp 



z — uz 

6 4 



dz . 



(C.3) 



(C.4) 



and C is any contour that starts at the point infinity with argument arg 2; = — | and ends 
at the conjugated point. Notice that we have used slightly different variables than in [20] . 
In order to prove the relation (IC.ip . we need to show that for any t G M, 



W'it) = - 



/oo 
P(2^2/3^^ 2-^/H)dm . 
-00 



(C.5) 



From a quick examination of (II. 4p and (IC.2p . one sees that (IC.SP holds given the following 
lemma, which is subsequently proven. 
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Lemma C.l. For any Xi, X2 > 0, m G M, and t G M, 
di 

Proof. Consider first the right hand side of (1C.6I) . Since 



(m, t) = 2-^/3^(21/^x1; -2-'/h, 2-2/3^)^(21/3x2; 2-'/h, 2-'/^m). 



(C.6) 



Ai{z) = — I e3*'-*^dt, zeC. 
2'Ki 



we have 

^(2i/3x;-2-^/3t,2-2/3m) 



e 2 



which after the change of variables C, = 2 ^^^r] — 2 '^^^t, simphfies to 

2-2/3 ,3^„ 



7/;(2i/3x;-2-^/3t,2-2/3m) 



-6 24 ' 4 / r/e 



The right hand side of (IC.6P can now be written as 

2-^V3^(2i/3xi; -2-^/3t, 2-2/3^)^(2^/3x2; 2-^/3^, 2-2/3m) 




(C.7) 



c JC 



For the left hand side of (lC.4p we have 

m Tfi 
^{x2 + - + y] t)^'(xi + y + 



1 



47r2 




J^-Cl(^2 + f)-C2(xi + f)-|(Cf-C|)-S/(Cl+C2)^^^^^, 



and 



Tfi Tfl 



1 



47^2 



Plugging these formulas into (IC.3P and interchanging differentiation and integration, we 
obtain 



$X2Xi("Z,t) 



^vr2 .,0 




e 6 



-Cl(x2 + f)-C2(xi + f )-|(C?-C|) 



£ 7£ 



X e-^(^^+^^) (Ci - C2) + e^^^^"^^^ (Ci + C2) dCidC2 



dy. (Ci 



The function 

/(y):=e-^K^+^^)(Ci-C2)+e^(^-'^^)(Ci + C2), Q e C, (17^(2, 
is integrable over [0, 00) since for G £ we have Re > 0, and if Re Ci < Re C2, then 

\f{y)\<CRe{C2)e-y'''^^'-^^\ 
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On the other hand for Re Ci > Re (2 we have for 5 > 

<CRe(Ci)(e-2^^^(^^) + e-^^), 
hence integrabihty follows. Since for Ci 7^ C2; 

Ci - C2 Ci + C2 4C1C2 



f{y)dy 



C1 + C2 C1-C2 cf-cr 

we can use Fubini's theorem in (IC.SP to obtain 







27r2 Jc Jc Ct - Ci 

This integral appears singular, but one of the contours C could be deformed from the begin- 
ning such that Ci and (2 never coincide. Now differentiating with respect to t (and exchanging 
integration and differentiation) gives 

|$.,.,(m,-t) = ^ |^|^CiC2e^-^^(--^f)-^^(-+f)+i(^?-^^^)rfCiC^C2. 
Comparing the last line with flC.7l) . we have shown 

|-$,,,,(m, -t) = 2-5/3^(21/3:^1; -2^^/% 2~''l''m)^{2}I^X2\ 2~'"^m). 
The change of variables t ^ —t gives (1C.6I) . □ 
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